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1 Introduction 



A well known [1] problem in high-temperature (T) field theory is the breakdown 
of the conventional perturbative expansion at some order in the coupling constant 
(g). This happens because in the regime T gT ^> m , where mo represents any 
intrinsic zero-temperature masses in the theory, the relevant cutoff for infrared (IR) 
singularities in loop diagrams is the thermal mass (~ gT) rather than m . Higher loop 
diagrams then accumulate powers of g in the denominator which can compensate for 
the usual factors of g in the numerator coming from the Feynman rules. Therefore, 
to compute consistently to a given order in g, we have to take into account all the 
relevant higher loop graphs — these usually form an infinite set. 

A practical solution is to resum the perturbation series by systematically includ- 
ing [2,3] all lower order radiative corrections that are significant (like the thermal 
mass) in higher order calculations. For gauge theories, the required resummation 
of the perturbative expansion into an effective expansion was developed recently by 
Braaten and Pisarski [3] to compute the gluon damping rate to leading ( ~ g 2 T) 
order. Subsequently, the effective expansion has been used to compute many other 
quantities [4] . In all these applications , only one-loop diagrams in the effective theory 
were considered. 

To go beyond leading order, one must compute two-loop (and higher) diagrams 
in the effective expansion. Since this is a tedious exercise in gauge theories, I will in 
this paper deal with a toy model — the g 2 4> 4 theory — in order to explore some of the 
technical aspects of higher loop calculations within the resummation program. As 
will be discussed in Sec. 2, for this model only the self-energy has to be resummed 
while the vertex can still be treated perturbatively as in the bare theory [3]. A two- 
loop calculation in the same model with partial resummation has been considered by 
Altherr [5] . More recently, a modified perturbation expansion for the model was pro- 
posed by Banerjee and Mallik [6] to enable the systematic calculation of the effective 
mass to higher orders. In [6] a mass parameter was introduced in the beginning and 
later determined by consistency conditions. 

The main difference between [5,6] and this paper is that here the imaginary-time 
formulation (ITF) will be used to perform the calculations whereas the real-time 
formulation (RTF) was employed in [5] and [6]. In the ITF the diagrammatics is 
the same as at T = and the power-counting of IR divergences is extremely simple. 
These conviniences of the ITF will be exploited to give a careful account of all the 
diagrams that can contribute to a given order in g towards the self-energy. Also, 
instead of introducing a mass parameter as in [6], the resummation will be done in 
stages so as to make it easier to identify the relevant diagrams and ranges of momenta 
which can contribute to a particular order in the coupling constant. As an example 
of an explicit calculation, I will determine the thermal mass and damping rate, for 
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bosons at zero momentum, up to order g 3 . The results will be compared with those 
obtained in the RTF. Of course, as the scalar model is quite popular, some of the 
formulae and results obtained in this paper, especially in Sec. 2, may be found in other 
publications [3, 5-8 ]. 

The plan for the rest of the paper is as follows : in Sec. 2 I will set up the notation 
and perform the first stage in the resummation of self-energy diagrams. This includes 
only one-loop one-particle irreducible (1PI) diagrams. At this stage, the self-energy is 
momemtum independent so the induced thermal mass is easily obtained to order g 2 . 
Although individual higher-loop 1PI self-energy diagrams in the effective expansion 
seem to contribute to this order, it is demonstrated that their sum does not. Thus 
perturbative computability is maintained in the effective expansion. The one-loop 
4-point function is also considered in order to explain why the vertex corrections can 
be treated perturbatively. In Sec. 3, the effective lagrangian of the previous section 
is used to perform the next stage of the resummation, which includes both one and 
two- loop diagrams. The thermal mass is obtained up to order g 3 . Again, the sum 
of higher-loop diagrams is shown to cancel at this order. The imaginary part of the 
self-energy is also computed to determine the damping rate to lowest order. The 
conclusion and a summary is in Sec. 4, while the appendix contains some technical 
details. 

2 One-loop 

The starting point is the following lagrangian for a hot scalar field (i.e. the intrinsic 
mass has been set to zero) 



The lagrangian has been written in D dimensional Euclidean space, where D = 
4 — 2e and /i is the mass parameter of dimensional regularization. The renormalization 
counterterms, which have not been displayed, will be determined in the minimal 
subtraction scheme [9]. In the imaginary time formulation of finite-temperature field 
theory [10,11], the information about the temperature (T) is contained in the energies 
which are now discrete; for bosons, p° = 2irjT , where j is an integer. The only change 
from the zero-temperature Feynman rules is then in the replacement (much of the 
notation is similar to [12] ) 



£ = - 




(2.1) 




(2.2) 
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The sum over the discrete frequencies inside loops is most efficiently performed 
using the "Saclay" method [12]. Real-time amplitudes are then obtained by analyti- 
cally continuing [10] the external energies, p° — > —iou. Let A(K) represent a bosonic 
propagator with mass M and momentum K 2 = (k ) 2 + k 2 , 

For the lagrangian £ , the massless propagator will be denoted as A (K) = 1/K 2 . 
In the Saclay method, the propagators inside loops are replaced by their spectral 
representations 



A(K) = / dr e A(r, k) 
Jo 

A(r,k) = (l/2E k )[(l + n k )e- EkT + n k e EkT ] (2.4) 



Here, El = k 2 + M 2 and n k = l/(exp(E k /T) - 1) is the Bose-Einstein distribution 
function. The expression (2.4) for the noncovariant propagator is valid for < r < 
1/T and is defined to be periodic in r with period \/T outside that range. By using 
the spectral representation of the propagators, it is trivial to do the frequency sums 
followed by the r integrals, leaving only the integrals over spatial momenta to be 
performed [12]. 

The main calculations in this paper will focus on obtaining consistently the pole 
of the effective propagator, l/(-P 2 — H(p°,p)), where P M = {p°,p) is the external 
4- momentum and II is the 1PI self-energy. For the theory described by (2.1), the 
diagram in Fig. la can now be evaluated as described above to determine the self- 
energy to lowest order 



W) = -^-Tr k A (K) 

g 2 ^ ( r d D ~ l k 1 r d^k n k \ 

2 \J (2tt) d - 1 2k J (2tt) d - 1 2k)' 1 ' 
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The first integral in (2.5) is the self-energy at T = 0. It vanishes in dimensional 
regularization [9], so there are no ultraviolet (UV) divergences to this order. The 
second integral in (2.5) represents the matter contribution and is UV finite because 
of the Bose-Einstein (BE) factor. Putting e = then gives the result 

q 2 T 2 

MP") = -^J-- (2-6) 

In this paper, the induced thermal mass, m, is defined as the real part of the 
pole of the Minkowski propagator at zero momentum (p — 0). Since the self-energy 
to this order is independent of momentum, one gets 



m 2 ee -n = 9 —. (2.7) 



To systematically include the effects of this thermal mass, the term |m 2 2 is 
added and subtracted [2,3] from (2.1) to define a new effective lagrangian 

C 2 = (Co + \m 2 <f) - \m 2 <f. (2.8) 

The subscript '2' on L 2 is used to remind us that the new lagrangian now describes 
a theory with tree level mass m with (m/T) 2 ~ g 2 . In (2.8), the quantity in brackets 
defines a lagrangian with free propagator A 2 (K) = l/(K 2 + m 2 ). The subtracted 
term is treated as a new 2-point interaction (Fig.lc) of order g 2 . The shifting of 
terms in C to form C 2 corresponds to a resummation of the perturbative expansion. 

The next step [3] is to use the effective lagrangian (2.8) to recalculate the self- 
energy. In addition to Fig. la, a contribution from the new vertex (Fig.lc) must also 
be included, 



,2,,2e 



B 3 (^) = m 2 -^TT k A 2 (K) 

2 g 2 m 2 fAna 2 Y , . g 2 a 2t f dP^k 2n k 



2 (4tt) 2 \ m 2 J v ' 2 J (2tt) d - 1 2E k 



where now E% = k 2 + m 2 . The second term in (2.9) is divergent as e — > 0. This 
UV divergence is similar to that in T = field theory. The only difference is that 
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as a consequence of the resummation the thermal mass has been introduced into the 
perturbative calculations. This makes the above divergence temperature dependent 
[6], albeit in a trivial way — the structure of the divergence is the same as at T = , 
with the intrinsic mass itlq replaced by the thermal mass m. Therefore the structure 
of the mass counterterm will still be the same as for T = 0, ensuring that the theory 
is renormalisable even though the counterterms are temperature dependent. (See 
however the discussion following (2.17)). Expanding the divergent term near e = 
gives 



— 7 — -T7- + finite terms $(g 2 m 2 In sL.) (2.10) 

2 (An) 2 e vy m 2> v 1 

where, since m 2 ~ (gT) 2 , the finite term is $(g 4 lng). The mass counterterm ver- 
tex (Fig. lb) is thereby fixed to be —g 2 m 2 /327i 2 e at lowest order. The one-loop 
renormalised self-energy in the theory defined by C 2 is then (See A. 4) 



3 4tt 2 Jo E k 



= 3m 3 /irT + $(g 4 \ng). (2.11) 



The corrected thermal mass, M 3 , is given by, 



M 2 = m 2 - W 3 en = m 2 (l - 3m/nT) + $(g 4 In g). (2.12) 

Calculating Fig. la using the propagator A 2 (K) is equivalent to summing the 
infinite set of "daisy" [13] diagrams of Fig. 2 evaluated with the massless propagator 
Aq(K). This interpretation follows once A 2 (K) = l/(K 2 +m 2 ) is expanded in a Taylor 
series in m 2 about m 2 = 0. Each of the diagrams of Fig.2 (for iV > 1) is infrared 
divergent in the theory Co but their sum is, as we have seen, infrared finite. Thus 
summing an infinite set of IR divergent diagrams has given an IR finite correction 
of order g 3 to the mass. The nonanalytic (in g 2 ) behaviour of this correction is a 
sign of its nonperturbative nature (infinite resummation) when viewed in terms of 
the original lagrangian (2.1) [5,11]. 

Are there any other diagrams in the effective theory (2.1) which can contribute 
terms of order g 3 to the self-energy ? The answer, at first sight, is yes. Even in the 
effective theory, there are infinitely many diagrams, other than those in Fig.l, which 
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can contribute at order g 3 but, fortunately for the consistency of the resummation, 
their sum is of order g A or higher! Let me term such diagrams "irrelevant" since 
eventually their finite contributions to the present order in g cancels , though they 
might be relevant for the UV renormalisation of the theory. 

Examples of irrelevant diagrams at order g 3 are given in Fig. 3. Each of the 
diagrams there is ${g 2 ) 2 (^) ~ g 3 - The factor (g 2 ) 2 comes from the vertices while 
the \ jg comes from the bottom loops as their IR singularity is cutoff by the thermal 
mass. The simplest way to deduce the factor of 1/g is to note that the IR behaviour 
of bosonic propagators in loops is dominated by the j — term in the frequency sum 
(2.2). That is , to get the leading IR behaviour of a diagram, set all the internal 
energies to zero and then take the m — > limit. Consider for example Fig. 3a. Its IR 
behaviour is 



d 3 k f d 3 q 



g*Tr k Tr q (A 2 (K)[A 2 (Q)] 2 ) ~ 9* J „ J 



k 2 + m 2 J (q 2 + m 2 ) 2 
% 4 )(l)(l/<?) = <% 3 ). (2.13) 



However, the sum of graphs in Fig. 3 ,with the proper combinatorial factors, is (using 
eqn.(2.11)) 



^Tr,[A 2 (g)] 2 f^Tr fe A 2 W - m 2 ) = (n 3 ) Tr q [A 2 (Q)] 2 



' % 2 )( 5 3 )(l/<?) = i% 4 ).(2.14) 



The UV divergent parts of course cancel only when all the relevant two-loop and 
counterterm diagrams are summed ( Sec. 3). Similarly, although each of the daisy-like 
diagrams shown in Fig. 4 is $(g 2 ) 3 (l/ g 3 ) ~ g 3 , their sum is easily shown to be 



-^Tr,[A 2 (g)] 3 ^Tr fe A 2 (X) - m 2 J ~ % 2 )(l//)(^) 2 = % 5 ). 

(2.15) 

In general each of the daisy-like diagrams in Fig. 5 with a fixed number iV > 2 of 
"bubbles" (Fig.la) + "blobs" (Fig. lc) is ?% 3 ) but their sum is 
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-gii Tr g [A 2 (Q)j ^ — ^— - 



p=0 

Tr, [A 2 (g)]^ L 2 - 9 -^f Tr k A 2 (K) 



N 



~<% 2 )(i/^-%^) = <V +3 ). 



N+3\ 



(2.16) 



Thus the set of all daisy-like diagrams with iV > 2 is completely irrelevant for the 
calculations in this paper which will be performed up to order g 4 . It is clear from 
the above analysis that the presence of the 2-point interaction (Fig.lc) is essential. 
Recall that it was introduced (2.8) to keep us in the same fundamental theory while 
performing the resummation. We see now how, in the cancellation of contributions 
from the infinite set of daisy diagrams, it prevents an overcounting of diagrams. 

The N — 1 daisies of Fig. 3 which seem to be relevant at order g 4 will be discussed 
further in the next section. It is left as an exercise for the interested reader to verify, 
using the simple power counting rules for IR singularities illustrated above, that any 
other 1PI self-energy diagram is individually of order g A or higher. 

To summarise, the thermal mass-squared including all subleading corrections of 
order g 3 is completely given by (2.12). 

So far, all the results have been written in terms of the renormalised coupling 
g. The 'physical' coupling is determined by evaluating the diagrams of Fig. 6 on shell 
, which corresponds to soft (~ gT) external momenta. Using the by now familiar 
power counting, it is seen that each of the diagrams is $(g 3 ) and hence the radiative 
correction to the basic 4-point vertex is down by a factor of g [7]. Therefore, vertex 
corrections can be treated perturbatively instead of resumming the corrections to 
form effective 4-point vertices. Contrast this with the thermal mass , m, which is of 
the same order as the bare inverse propagator Aq 1 (K) for soft momenta and therefore 
has to be resummed. In the language of [3], for the scalar theory, the only " hard 
thermal loops" are in the self-energy. In this paper, all the results will be left in terms 
of the renormalised coupling g. 

To obtain the complete effective lagrangian to order g 3 the vertex renormalisation 
counterterm is needed. This is determined as usual by calculating Fig. (6a) ( plus the 
usual crossed diagrams ) at T=0. Including first all the counterterms in (2.8) gives 
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L2 t „2^2 



C 2 ^C 2 = (Co + m 2 2 /2) + ^ 
- ( m 2 2 /2 + 



9 ™ M + 9"» i 



2..2e / 3^2 l 



2! V(4vr) 2 2e 
2 g 2 m 2 1 \ 



4! r V( 47r ) 22e , 



(2.17) 



2! (4vr) 2 2ey ' 



Included in £' 2 is the counterterm (Fig. 7a) for loop corrections to the 2-point 
interaction. Note that the net lagrangian does not contain temperature dependent 
counterterms though pieces of it do because of the resummation [16]. 

Then, as before, the effects of the thermal mass to order g 3 are included by 
shifting the mass term in (2.17), 



C 3 = (C + M 3 V/2) 



M 2 072 + 




^2 „2 



2! (4tt) 2 2e, 



+ 



9 2 ^ 

IT' 



% 2 r 



(2.18) 



Notice that for consistency, the mass in the UV counterterms has also been shifted 
to M 3 in order to cancel the divergences in loop calculations. The lagrangian (2.18) 
will be used in the next section to obtain the self-energy up to order g 4 . Strictly 
speaking, the further resummation to obtain £ 3 is unnecessary as the $(g 3 ) correction 
to m 2 is a perturbative correction, just like the $(g 3 ) correction to the 4-point vertex. 
However, no harm is done by this additional resummation, all that happens is a 
redistribution among the diagrams of the next correction at order g 4 , as we will soon 
see. 



3 Two- loop 

The basic diagrams that must be considered to evaluate the self-energy to order g 4 
using £ 3 are in Figs. 1,3,7 and 8. Before delving into the calculations, let us make 
some observations. Just as for T = 0, the sum of graphs in Fig.l must be UV finite. 
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Also, as for T = 0, Figs. 7b and 7c are the mass and vertex counterterm diagrams that 
are required to cancel the sub divergences arising from the two-loop diagrams Fig. 3a 
and Fig.8 . A 'new' [6] feature of the resummation is the counterterm of Fig. 7a which 
is needed to cancel UV divergences generated by loop corrections (Fig. 3b) to the 
2-point vertex (Fig.lc). 

The sum of diagrams in Fig.l, evaluated from the lagrangian £ 3 is, 



rr?V) 



Ml 



n 2 n 2e a 2 M 2 1 



2 (4tt) 2 e 



= Ml-m 2 



3m 



m 



m\ 2 3 f m\ 2 



'-^-^{t) Ht) --At) { c >-T« 



+ 



1 /gm\' 

2 V4W 



]nf^U(l-7-) 



+ % 5 lny), 



(3.1) 



where ^ E is the Euler constant and C\ = \{2^ E — 2 ln47r — 1). Eqns.(2.12) and (A. 4) 
were used to get the final form (3.1). 

Fig. 3 contributes 

n?V) = ^ Tr, [A 3 (Q)f Tr fe A 3 (K) - M 3 2 ) . (3.2) 

Since (^ 2 /2)Tr fe A 3 (^) = M 2 + ${g A In g), therefore the finite part of (3.2) is?% 5 ln#). 
In the last section, working with C 2 , it was shown that the same diagrams sum 
to ^((? 4 ). The sum has now been pushed to higher order because of the further 
resummation performed to obtain £ 3 . This is an example of the 'redistribution' 
mentioned at the end of the last section — the 'lost' contribution from Fig. 3 has 
been picked up by Fig. la: this is indicated in (3.1) by the presence of the factor 
(Ci — 3/27r) rather than C±, the latter factor being the contribution if £2 were used 
in calculating Fig.l. The diagrams of Fig.3 are thus only needed to complete the UV 
renormalisation of the theory which, as usual, is performed loop-wise. 

The only other graph that is relevant for discussion is given in Fig.8; it will be 
considered later. The daisy-like diagrams (with N > 2) were already shown in the 
last section to be irrelevant at ^(g 4 ). However now a new infinite set of graphs must 
be analysed, the simplest of which are shown in Fig.9. Each of the diagrams in Fig.9 
is of order g 4 by power counting but their sum is clearly $(g 6 Ing). Extending diagram 
(9a) by adding a bubble (or blob) to its top gives a graph of order g 5 . So one only 
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needs to consider adding N number of bubbles + blobs to the middle loop and M 
bubbles + blobs to the bottom loop of Fig. 9a to create the general " cactus" diagram 
of order g A shown in Fig. 10. It is sufficient to show that the sum of all such cactus 
diagrams is of order g 5 or higher : first consider Fig. 10 with the bottom loop and 
its M attachments fixed in a particular configuration. Then any subdiagram (with 
fixed N > 0) above the bottom loop is precisely a daisy diagram and these have been 
shown to sum to $(g 5 ) at most. For any M > the bottom loop in Fig. 10 contributes 
a factor g 2 (l/g) = g. Hence the sum of all possible cactus diagrams is at most of 
order g 6 . This completes the proof. 

Adding bubbles or blobs to Fig. 8 creates diagrams like those shown in Fig. 11. 
These sum to order g 6 lng. All other diagrams are individually of order g 5 lng or 
higher. 

Having accounted for all the relevant diagrams, let us return to some explicit 
results. The counterterms in Fig. 7 contribute 



n 



(7) 



v i- 



<?V e 



2e 
3^ 2 



+ 



(An) 2 2e 



gM 3 

An 

Tr k A 3 (K) 



2e 



Tr fc [A 3 (K)f 



(3.3) 



Summing (3.2) and (3.3) gives 



n 



(3,7) 



2e (An) 2 

_g^Ml 
2e (An) 



( lE - 1) - In 



'Ann 2 

7m[ 



3g 4 Mf 
4e 2 (An)' 



+ finite terms $(g 5 In g), 



(3.4) 



The first line in (3.4) is a nonrenormalisable temperature dependent infinity gen- 
erated by diagrams (3a) and (7c). It will cancel [8] when the two-loop overlapping 
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diagram of Fig. 8 is added. This last diagram is the only relevant diagram which 
depends on the external momentum (p°,p), 



nfV,$ = ^Tr k Tr q [A 3 (K)A 3 (Q)A 3 (P-K-Q)} 

= G (p ,p) + G 1 (p ,p} + G 2 (p ,p) (3.5) 



where 



Go(p°,P) = Jd[k,q]S(E k ,E q ,E r ) (3.6) 
Gi(p°,p) = 3jd[k,q]n k [S(E k ,E g ,E r ) + S(-E k ,E q ,E r )] (3.7) 



G 2 (p°,P) = 3 / d[k, q] n k n q [S(E k , E q , E r ) + S{-E k , E q , E r ) + S(E k , -E q , E r ) - S(E k , E q , 

(3.8) 

with the definitions 

S(E k , E q , E r ) = 



Jp° + E k + E q + E r -ip° + E k + E q + E r§ 



iri . g 4 fi 4e d D ~ l k d D ~ l q 1 
d[k,q\ = 



3! (2tt) d - 1 (2tt) d - 1 8E k E q E r 
r = \k + q — p\ 

">2 ;2 i i\/r2 



Ef = l 2 + Mi , l = k,q,r 



and ni the usual BE factor. The real-time retarded self-energy follows by making the 
analytic continuation p° — > — iuo + £ with £ = + [10]. Then the prescription 
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gives the real and imaginary parts of the diagram [14]. Consider first the real part. 
Since G does not contain any Bose-Einstein factors, it must be the expression for 
diagram (8) obtained using T = Feynman rules and with the energy integrals done. 
In covariant (i.e. with P 2 = —uo 2 + (p) 2 ) notation one gets [15] 



ReG (P 2 ) = Re 



4 4e 

9 K 
3! 



d D k r d D q 1 1 1 

[2n) D l (2tt) d K 2 + Ml Q 2 + Ml (K + Q - P) 2 + M| 



-9 A M. 



1 3-2te 2 
-t; H h - In 



4 (4vr) 4 

+ finite terms $(g 4 M 2 )(P 2 /M 2 ) 2 



4lTfi 2 



4 (4tt) 4 6e 



(3.10) 



For soft external momenta (P 2 ~ m 2 ), the region of interest, the finite terms are of 
order g 6 and so do not contribute to the self-energy at order g 4 . 

Gi represents the mixing of the T = piece from one loop with the T^O piece 
from the second loop. This is clear from the expression (3.7) which contains only one 
BE factor, making one of the loop integrals UV finite while the other loop integral has 
a UV divergence. Specialising to the case p = in order to do the angular integrals, 
gives (See Appendix) 



Re Gi(-iw,0) 
where F 



and Foiui 2 ) 



F + F 1 + F 2 {uo 2 ) 

9* I e p(M 3 ) 1 
2 (An) 2 e 

g 4 If°(M 3 ) ( 4^ 

T^4^ !!U - 



Ml 



JO -C/fc JO tin I 



8(2tt) z 



3 
13 



X 4 



4k 



(3.11) 
(3.12) 



(3.13) 
(3.14) 



withX± = cu 2 - (E k + E q + E k±q f cu 2 - (E q - E k + E k±q ) 



The temperature dependent infinity F is actually independent of the external mo- 
menta p^ 1 and cancels precisely against a similar term found earlier in eqn.(3.4). 

Finally, G 2 contains a BE factor for each loop and so is UV finite. It however has 
a logarithmic IR divergence asm,u^ 0. One obtains 



H(u 2 ) = ReG 2 (-iuj,0) 



8(2tt) 4 Jo 



00 „ kn k 
ak- 



E k Jo 



(3.15) 



where 



Y. 



± 



(E k + E q + E k±q ) 



LO 



(E q — E k + E k±q ) 



X 



^ 2 — (E k — E q + E k ± q ) 2 J 2 — (E k + E q — E, 



(3.16) 



The sum of all UV divergent terms from eqns.(3.4, 3.10, 3.12) gives 
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2 \ 2 



g 2 \ 2 P 2 



[16^ J M "{^-Ae)-{Tb) 2le- ^ 
These are cancelled by the two-loop wave-function and mass renormalisation coun- 
terterms 



\ 16?r 2 y 24e 



+ 



Mf 



2 \ 2 / 1 
9 W 1 



16tt 2 



1 



2e 2 4e 



(3.18) 



The temperature dependent UV mass counterterm in (3.18) is i)(g 6 ). As before, it 
will precisely compensate [16] the i)(g 6 ) temperature dependent UV counterterm for 
the 2-point interaction (the diagrams that require the latter counterterm are formed 
by adding a blob to Figs. 3a , 7b, 7c and 8). 

The real part of the renormalised self-energy for £ 3 is therefore (3.1,13-15) 



R(uo 2 ) = Re Yl r 4 en (-iuj, 0) = IT^ + F 1 + F 2 (u 2 ) + H(uo 2 ) + terms of order (g 5 In g). 

(3.19) 
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This expression contains all corrections at order g 4 . It also contains some effects at 
order g 5 and higher in the energy dependent terms F 2 and H . For general uj, the 
expressions F 2 (ou) and H(oj) are too complicated to evaluate in closed form. However 
since only contributions to $(g 4 ) are required, something can be said. Note that 
because of the explicit factor of g 4 , it is only necessary to identify the IR behaviour 
of the integrals in eqns.(3.14-15) to obtain information about the leading contibution. 
Clearly, the uj dependence of any IR behaviour in F 2 or H can only be possible for uj 
soft (~ m). Consider first F 2 (oj 2 ). It is easy to see that the logarithmic IR singularity 
as m, uj — > is due to the second factor in the term X_ in the region k > q. So one 
is led to investigate the piece 



00 xdx f x ydy 

fir / — — In 



E x Jo E y 



{E y + E x - V - E x f - a 2 , (3.20) 



where I have factored out T 2 and defined the set of dimensionless variables {x, y, a, a} 
by scaling the quantities {k, q, m, uj} respectively by 1 jT . Now, for a — > the estimate 
(E y + E x _ y — E x ) 2 ~ d(a A )(l/y + 1 /{x — y) — 1/x) 2 holds. From this it can be deduced 
that for a ~ a n ~ g n the 'leading log' contribution from (3.20) goes like 2 In a 2 for 
n > 2 and like nlna 2 for n < 2. For H(uj 2 ), the logarithmic IR singularity is caused 
by the extra BE factor while the magnitude of its contribution is controlled by the 
m(y + /y_) term. Writing 



^ = -(8xy) 2 (x ± yfa 2 + (a 2 - a 2 ) 2 Z ± (3.21) 



where Z± = 16 (x ± y) 2 (x 2 + y 2 ) + (xy) 2 + a 2 (x 2 + y 2 ± xy) 
+8(a 2 - a 2 )(x 2 + y 2 + a 2 ± xy) 
+(a - a ) , 



shows that the first term in (3.21) dominates for a ~ a, while the transition to more 
complicated behaviour is again at a ~ ga ~ g 2 . 

To get the full g 4 dependence from F 2 and H, the constant under the leading 
g 4 In g contribution is also needed. This can be done for specific values of uj when 
it is possible to isolate clearly the q 9(g 4 ) pieces from the partial higher order effects. 
A calculation, which is sketched in the appendix, gives (with corrections at $(g 5 ) 
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omitted) 



F 2 (0) 
H(Q) 

F 2 {Ml) 



A 



7T 



A 



A 



rn 



H{Ml) = A 



hi ( - ) +3.48871... 



4 ln (?) + - 54597 - 
^ln(^) 2 + 4.52097... 



(3.22) 



from which follows 



F 2 (0)+#(0) = A 



F 2 {Ml)+H{Ml) = A 



m 



lnf — ) + 5.3025. 



m 



hi ( - ) + 5.0669... 



(3.23) 



where A = — (gin/An) 2 . Surprisingly, from (3.23 ) it appears that the coefficient in 
front of the total q 9{g i In (?) contribution to the real self-energy (3.19) from the energy- 
dependent part is the same on-shell as for zero external 4-momentum. Let me now 
proceed to the determination of the pole of the effective propagator. The complex 
pole, Q, at zero-momentum (p = 0) is the zero of the equation 



- Q 2 + Ml - U r 4 en (-icu, 0) = 0. 



(3.24) 



Since Im n 4 (— ifl, 0) is ^(g 4 ) (see later), then by writing Q = uj — ry, the real part 
of the pole up to $(g A ) is determined by 



uj 2 + Ml - R{uj 2 ) = 0. 



(3.25) 



The above equation may be solved by iteration (See Appendix) to give the thermal 
mass M 4 up to order g 4 , 
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M 4 2 = Ml- R{Ml), (3.26) 

where the right hand side must be expanded up to order g 4 . Using the values for 
F 2 (M|) and #(Mf) given in (3.23), together with eqns. (3.1), (3.13) and (3.19), in 
(3.26) above gives the final answer 



Ml = m 2 ( 1 



3m\ f gm\ 2 
ttT) + V 4W 



I In 

2 



T 2 , mX 2 

+ 2 In — 

47T/i 2 VT 



/(/mV 



(3.27) 



where a = 14.1416 and m 2 is defined by eqn.(2.7). Taking the square root of 

the above expression gives the complete thermal mass up to order g 3 . 

The imaginary part [14] of the self-energy to order g 4 is due only to the two-loop 
diagram of Fig. 8. From (3.5-9), it is relatively simple to obtain the imaginary part at 
zero-momentum and on-shell , 



im ni 8) ( 



-iM 3 ,0) 



16 \2tJ 



g 2 m 2 
32tt 



Jo E k Jo E q 



qdq 



+ ^ 5 )- 



(3.28) 



The result, as expected on general grounds [12], is positive. The imaginary part could 
also have been obtained directly, without using the prescription (3.9), by keeping the 
full logarithm in (3.14-15) instead of only its principal value. Finally, the damping 
rate follows from (3.24), 



7 



Im nfHM^o) 



g 2 m 
64tT 



2M 3 
+ #(g 4 



(3.29) 
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4 Conclusion 



To summarise, the effective expansion created by a resummation in the original the- 
ory was used to obtain the zero momentum pole of the propagator in the energy 
plane consistently to order g 3 . Working order by order in the effective expansion, 
the relevant diagrams were identified and perturbative computability was shown to 
hold by explicitly verifying the cancellation of contributions from an infinite class of 
diagrams. 

The first resummation of self-energy diagrams to get the effective lagrangian £ 2 
was essential because the thermal mass at lowest order, to, is as large as the inverse 
massless propagator at soft momenta. That is, the thermal mass could not be treated 
as a perturbation. As pointed out in the text, the second resummation to form C 3 
was not really necessary since the order g 3 correction to to 2 is a perturbative effect. 
The consequence of the second resummation was simply to change the individual con- 
tributions from some of the diagrams at "$(g 4 ). In particular, whereas the diagrams in 
Fig. 3 would have been relevant if we had continued using £ 2 , they became irrelevant 
when £ 3 was used — this 'lost' contribution was compensated by new subleading con- 
tributions from the diagrams in Fig.l. In short, though the reader could have been 
spared any mention of £3, nevertheless the author feels that some insight into the 
effects of resummation was gained by the exercise. 

Let me now make some comparisons with results in the literature. In [6], a mass 
parameter was introduced in the beginning and the effective mass was defined by 
requiring that the corrections to the free inverse propagator vanish at zero external 4- 
momenta. Translating that into the language of this paper, simply amounts to using 
R(0) on the right-hand side of (3.26) rather than i?(M|). In general, by definition, 
this does not give the pole in the propagator. From (3.23) we see that the difference 
between the two definitions shows up in the constant under the logarithm at order 
g 4 lng (in the appendix I explain the difference). I do not know however whether it is 
a coincidence that the coefficients of the 'logs' in (3.23) are the same. The imaginary 
part of the self-energy of course vanishes for u = 0, as is apparent from (3.5-9) or 
from more general arguments [12], but is given on-shell by (3.28). 

Clearly, any quantity calculated in either the RTF or ITF must give the same 
result even though some of the intermediate expressions may look different because 
of the differences in approach. For completeness, I have checked (using the real-time 
expressions found in [5] and [6]) up to order g A that the two formalisms give identical 
answers for the pole of the propagator; and also when used to calculate the effective 
mass as defined in [6]. 

For gauge theories, both the one-loop self-energy and vertices must be resummed 
into effective quantities [3]. Since these quantities are momentum dependent, the 
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effective expansion is quite involved even at one-loop order. Nevertheless, one expects 
that some of the features of two-loop calculations studied here in a simpler context 
will also manifest themselves in gauge theories. 
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A Appendix 



(1) The basic expression appearing in one-loop diagrams is 

fi 2e TrA k (K) = J (M) + J|(M) (A.l) 

with 

/„(«)= ^/^ s ^ T _ = _^j r(-l + e) (A.2) 
For (M/T) <C 1, we have the expansion [13], 



(A.4) 

The expression /x 2e Tr [A(K")] 2 can be obtained by differentiating (A. 1-3 ) with 
respect to M 2 . 

2. To get (3.11-14). 

Consider the q integrals in (3.7) . For p = the only nontrivial angular integral in 

— * — * 

(D — 1) dimensions is for the angle 9 between k and q. Choose k to define the polar 
axis, and first do the trivial angular integrals for the q- variables (see [9] , for example, 
for the correct measure for the integrals in D dimensions). Then one is left with the 
following integral (t = cos 9), 



dq !_ 
Jo " 

where 



r^q^L(k,q) (A.5) 

JO h a 



L(k,q) 



1 , 1 



d 



t 2 y at 



(in 



(E q + E r — E k y 



LO 



{E q + E r + E k ) 



(A.6) 



oo 



The simplest way to proceed is to subtract the UV divergent part of (A.5). As q 
' §t — *■ 2/c/g. Subtracting and adding this term at the appropriate place in (A.6), 
substituting everything back in (A.5), (3.7) and then doing the obvious simplification 
gives the result quoted in the text. 

3. Solving Eqn.(3.25). 

Since Mf ~ $(g 2 ), and R(uo 2 ) ~ ^(g 4 ), a consistent way to solve (3.25) is by iterating 
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the lowest order solution uj 2 = M 2 . The next iteration gives the result in the text 
(3.26) while further iteration will give a correction at i)(g 6 ). Essentially then, the 
pole is determined by the self-energy on mass-shell. Now consider a Taylor expansion 
of R{oj 2 ) about uj 2 = 0. Since R has a logarithmic IR singularity as u,m — > (see 
text), therefore 



{u 2 ) n d n R{u 2 



For oj soft (~ m), the Taylor expansion is not an expansion in g as each term is of 
order g 4 . So one should expect R(0) and R(oj 2 ) to differ by an amount $(g 4 ) (see 
(3.23)). The same argument explains why it would be difficult to obtain the full order 
g A contribution to the self-energy at soft non-zero external momentum (p) by doing 
a Taylor expansion about p = 0. 

4. To obtain (3.22). 

To extract the leading "lng+ constant" contribution from the integrals appearing in 
(3.14-15), the following procedure is adopted : identify the terms which will contribute 
the lnm singularities, isolate them and then set m = in the regular terms to get 
pieces of the ^(g 4 ) contribution. Next, simplify the potentially singular terms and 
keep repeating the above procedure until all the "log + constant " pieces have been 
explicitly obtained. Consider for example F 2 {uj 2 = Mf). The term within braces in 
(3.14) is simplified (for oj on-shell) and written (replacing M 3 with m in the expressions 
F 2 and H ignores a correction of order g 5 ) as follows 



2n 



gin 

—gin 
-4k 



' q + k \ ( m 2 + q(q + k) + E q E q+ k 
<l /J I " T 2 " " 
m 2 + q(q - k) + E q E k _ q 



+ 



J^2 



+ 



(A.8) 

(A.9) 
(A.10) 



As discussed in the text, there are no IR singulrities in the region q > k , so one may 
set m = in the above expressions in that range of integration. For the q < k sector, 
(A.9) and (A. 10) both contribute 'logs' while (A.8) gives a finite piece in the massless 
limit. For the piece in (A.10),the q integral is easily done explicity, then the lnm piece 
isolated and the finite terms determined. For (A.9), multiply the argument of the 
logarithm by (m 2 + q(q — k) — E q Ek- q ) in the numerator and denominator, simplify, 
isolate the 'log' piece and set m = in the rest to get the constants. Collecting all 
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the terms gives 



F 2 (M|) = 



g 1 



6 64tt 2 



2 y T' 



In 2 



4/ 

7T Z JO 



"° xhix 
ax 



+ $(g*\ng) (A.ll) 



The integral in the final answer (A. 11) is a pure number. It may be computed 
numerically if required. The concise result is given in the text (3.22). 

For a different example, consider 



H(Mi 



g 



00 „ kn k 
dk 



dq— In 

o E a 



k 



k 



q 



(A.12) 



64tt 4 Jo "" E k 

where the symmetry of the integrand under k <-> q interchange has been exploited 
to restrict the range of one of the integrals. As the IR singularity is now due to the 
extra BE factor rather than the explicit logarithm, the 'log + constant' pieces can be 
obtained by using an arbitrary soft cutoff A (I thank R. D. Pisarski for suggesting 
this technique) to divide the region of integration for the k variable. For k < A, 
the appropriate approximations can be made ( e.g. n k = 1/E k ) to simplify the 
integrations . In the limit a — > 0, one gets for k < A in (A.12), 



T 



^lnS + ^lnd)+ fdtln 
4 y T' 4 y T' Jo 



1 + t 


hit 


1-t 


(f-f3)J 



+ tf(m)+tf(A) (A.13) 



For the region k > A, one can put the mass to zero because there are no lnm 
singularities. Next isolate the In A factor (for A — > 0) by doing a subtraction of the 
leading IR part of one of the BE factors, to get finally for k > A in (A.12), 



T 



7r 2 , . A. r 1 , r°° , x 
— In - + / dt / dx-^ 
4 1 Jo Jo e x — 1 e 



xt 



xt' 



In 



1 + t 



+ tf(m)+tf(A) (A.14) 



In the limit A — > , the sum of (A.13) and (A.14) gives the final answer for (A.12), 
with neglected terms of order g 5 . The cancellation of the In A terms in the sum 
removes the ambiguity coming from the cutoff. 

Similar considerations as above give 
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g 4 t 2 vr 

6 64tt 2 



(A.15) 
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32tt 4 
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4 T 2 r 1 



32tt 4 
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dx 



Ve xt 
tx 3 



l)lnf 1 + t + ^ 

arT \ 1 — t + ^ 2 y 



(x 2 + l)(x 2 t 2 + l) 



In 



4x 2 + 



(i+t+t 2 ) 



4x 2 + 



(i-t+t 2 ) 



(A.16) 



Again, if necessary, the constant integrals appearing above can be done numeri- 
cally to give the result in (3.22). 
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latter fact is transparent in the approach used in [6]. 
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Figure Captions 



Fig-l : 

Fig. la is the one-loop self-energy diagram, also called the "bubble". 
Fig. lb shows the ultraviolet mass counterterm while Fig.lc is the finite 
2-point interaction ( "blob" ) counterterm induced by the resummation. 

Fig-2 : 

Fig. 2 is the "daisy" diagram with N > 1 attachment of bubbles. 
Fig-3 : 

The the set of N — 1 daisy-like diagrams; Fig. 3a is the N = 1 daisy 
while Fig. 3b has an insertion of the finite 2-point interaction (blob) 
into the one-loop self-energy (bubble) diagram. 

Fig-4 : 

The set of N = 2 daisy-like diagrams. 
Fig. 5: 

A general daisy-like diagram with N > 1 bubbles + blobs attached. 
Fig.6 : 

Fig. 6a is the one-loop correction to the 4-point vertex, the diagrams in 
the crossed channels are not shown. Fig. 6b is the UV vertex countert- 
erm. 

Fig-7 : 

Renormalisation counterterms for second order calculations. Fig. 7a is 
the counterterm for Fig. 3b. The mass and vertex counterterms are 
Figs. 7b and (7c) respectively. 

Fig-8 : 

Overlapping two-loop self-energy diagram. 
Fig.9 : 

Diagrams which are individually of order g A but sum to higher order. 
Fig.10 : 

A general "cactus" diagram. N > bubbles + blobs are attached to 
the middle loop while M > bubbles + blobs are attached to the 
bottom loop. 

Fig.ll : 

More complicated self-energy diagrams which are individually of order 
g A but sum to higher order. 
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